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ABSTRACT 

The problem of formation of generic structures in the Universe is addressed, whereby first 
the kinematics of inertial continua for coherent initial data is considered. The generalization 
to self— gravitating continua is outlined focused on the classification problem of singularities 
and metamorphoses arising in the density field. Self— gravity gives rise to an internal hierar- 
chy of structures, and, dropping the assumption of coherence, also to an external hierarchy 
of structures dependent on the initial power spectrum of fluctuations. 

I. Introduction 

One of the ambitious goals of modern cosmology is to understand the morphology of 
large-scale structure as it manifest itself in the distribution of galaxies, but hiding its 
origin and key-elements of understanding in the distribution of dark matter. A detailed 
analysis of evolution models for large-scale structure is necessary to uncover the full 
morphology of the, according to present knowledge, dominating dark matter constituent 
in the Universe. 

The formation and evolution of generic structures is another fascinating variation of 
this theme, supported today by a well-developed mathematical literature on Lagrange- 
singularities and metamorphoses (see: Arnol'd et al. 1985, 1988), and sought for in a 
variety of areas such as for example morphological studies in biology, only to mention 
one of the fruitful applications of catastrophe theory (see: Thom 1975). 

Understanding the global morphology of structures with the help of generic singu- 
larities that arise locally in the dark matter density field is also the focus of this talk: our 
aim is to identify the building bloclis of large-scale structure. The problem here is that 
the possibility of classifying generic structure elements into a finite set of morphological 
germs is thus far restricted to continua which move under inertia; structure formation in 
the Universe is, however, a result of self-gravity. We shall approach this problem by first 
developing the principal picture of the much simpler problem of structure formation in 
an inertial continuum, then generalizing the framework to self-gravitating continua along 
the same track of argumentation. Indeed, we shall see that self-gravitating structures can 
be understood to a large extent on the basis of structures that arise in inertial continua, 
the gravitational problem, however, reveals large qualitative differences both in the fun- 
damental equations and in the possible richness of morphological manifestations. While 
the motion of inertial continua is well-imderstood, the research area of self-gravitating 
motions is still a sparsely charted terrain with many challenging problems lying ahead. 



II. Structure formation in inertial continua 



II. 1. Lagrangian kinematics of inertial continua 

In the Lagrangian picture of continuum mechanics we consider the deformation field: 

x = f(X,0 , (1) 

where X are Lagrangian coordinates labelling fluid elements, x are positions of these 
elements in Eulerian space at the time t, and f is the trajectory of fluid elements for 
constant X (vectors are written bold in this paper). The velocity field of the fiuid is the 
tangential field to this family of curves, v := f ; the dot denotes the total (Lagrangian) 
time-derivative along the curves. 

Inertial motion means that the acceleration of fluid elements vanishes along the flow 
lines, g := f = 0. Hence, the Lagrangian evolution equations can be written as: 

f = , (2) 

which we immediately integrate to obtain the general solution to this problem: 

f = X + V(X) (t-to) . (3) 

The integration constants have been chosen to have the initial positions coinciding with 
the Lagrangian coordinates, f (X, to) =: X; the initial velocity fleld is the second inte- 
gration constant, v(X,to) =• V(X). 

The velocity of each fluid element is constant along the flow lines, but can differ for 
different elements, i.e., the velocity field is in general inhomogeneous. We shall assume 
throughout most of the following considerations that the initial velocity field is coherent 
on some specified spatial scale on which the structure formation process is studied. This 
assumption we shall drop at the end of this talk. 



II. 2. Generic structure formation in inertial continua 

Consider, e.g., a one-dimensional sinusoidal velocity distribution. Continuum elements 
with higher velocity will catch up those elements with lower velocity and will overtake 
them. At the instant of overtaking the velocity field in Eulerian space is vertical, after 
that it will become multi-valued, i.e. three continuum elements arrive at the same 
Eulerian position. This situation is called three-stream Row. 




Figure 1: The one— dimensional motion is illustrated in the three— dimensional 'phase space' 
{X, X, t): trajectories of continuum elements (starting at X = X) trace the 'Lagrangian manifold' 
JO (defined by dx A dX = 0). In the case of inertial motion, there is a single folding of jO after 
shell— crossing; the projection of the critical points of the mapping X I — y X form the caustic in Eulerian 
space enclosing a three— stream flow. 

If the initial density field is ^o(X), then the density field in Lagrangian space at the 
time t changes inversely to the deformation of the continuum. In three dimensions this 
deformation is measured by the determinant of the deformation tensor (dfi/dXk), and 
the density in Eulerian space is: 

,M = ± ,{h^Mt) = ± ^rf)i''o^\^ r ' (4) 

^1 ^^\det[5ik + {t-to) ||^(h^[x,i],i)| 

where = hm[x, t] denote the m solutions of the inverse of the transformation (3); n is 
the number of streams. At the time of verticality of the velocity field, a singularity in the 



mapping (3) from Lagrangian to Eulerian space develops; images of these singularities lie 
on surfaces with infinite density in Eulerian space (the caustic). Caustics which enclose 
three-stream flows are called pancakes following Zel'dovich (1970) and Arnol'd et al. 
(1982). 

The total density in the pancake is the sum of the moduli of the individual densities 
of the streams according to (4). This is correct as long as the streams do not interact. 
The one-dimensional picture is illustrated in the space {x,X,t) in Fig.l. 

In two-dimensional continua the caustics associated with pancakes have the shape 
of sickles with cusped edges (compare Arnol'd et al. 1982). In order to understand 
the further evolution of pancakes, we consider the velocity fronts which are at each 
point orthogonal to the trajectory field. In Fig. 2 we illustrate what happens for an 
initially coherent and anisotropic velocity front. For further discussion of density field 
metamorphoses the reader may study the articles by Zel'dovich (1970, 1973), and in 
particular Arnol'd et al. (1982) for two-dimensional, and Arnol'd (1982), Buchert et al. 
(1994) for three-dimensional singularities and metamorphoses. 

The finite classification of structure elements associated with singularities of the 
mapping (3) is possible for space dimensions < 5 (Arnol'd 1982, Arnol'd et al. 1985, 
1988; see also the contribution of Shenglin Cao (1993) in the Proceedings of our last 
workshop). However, the mapping (3) defines a one-parameter (t) family of Lagrangian 
mappings only, if the motion is potential, i.e., the mapping (3) from Lagrangian to 
Eulerian space is a gradient mapping for each t. Arnol'd et al. (1982) and Buchert et al. 
(1994) give the local algebraic forms of singularities as determined by the eigenvalues of 
the deformation tensor and its derivatives. This way we are able to predict the formation 
and transformation of structure elements directly from initial conditions. 

The kinematics of an inertial continuum can be understood from the local degen- 
eracies of the deformation tensor: Following Zel'dovich &: Myshkis (1973) we write the 
density field (4) in a diagonal frame at a given point X for the single-stream potential 
motion, V =: — Vx<5, 



'^<^'''=(l-xc,(X))(l-x«X))(l-r7(X)) ■ ' 



where a > /3 > 7 denote the eigenvalues of the tensor (d^S /dXidX^). 

In the (non-generic) situation, where all the eigenvalues are equal (spherically sym- 
metric collapse), a continuum element degenerates into a single point. In the generic 
case, however, the largest eigenvalue a starts to dominate the collapse yielding a strongly 
anisotropic degeneracy of volume elements into surface elements orthogonal to the eigen- 
vector associated with a. After that the second eigenvalue P degenerates: the sheet-like 
element shrinks into a line element along the eigenvector associated with /3. Finally, 7 
degenerates, and the line element shrinks to a point-like matter concentration. 




Figure 2: The two— dimensional motion is illustrated in terms of velocity fronts. Initially, the 
anisotropic deviations of a coherent velocity front from an isotropic (spherically symmetric) front, 
which would collapse into a single point, is shown. Relative velocity vectors pointing inward move the 
fronts until they touch each other at a finite instant of time (a 'Legendre singularity' of the front forms 
corresponding to a 'Lagrange singularity' of the mapping X I — ^ X). After that time the fronts overlap 
enclosing a three— stream system of the flow corresponding to a, 'pancake' caustic in the density field. 
The pancake boundary breaks (shown here as the break of the velocity front) with edges connecting 
other regions of space: this is the origin of a global network structure. 

However, this picture is idealized, since after the sheet elements have formed, the 
mass shells cross freely and form a three-stream system with finite thickness. The prin- 
cipal picture of these subsequent degeneracies is illuminating: it explains locally the ex- 
istence of fiat sheet patterns, filamentary elements and matter concentrations in clumps. 
If the density field is coherent, this local picture can be extended into the surrounding 
space, since the eigendirections associated with the eigenvalues only vary slightly with 
position. Hence, a global picture arises with sheets, filaments and clumps: the pancake 
picture of structure formation (see: Shandarin &; Zel'dovich 1989). The density field 
around such a clump of matter is illustrated in Fig. 3. 



III. Structure formation in self-gravitating continua 



III.l. Lagrangian dynamics of self gravitating continua 

The Lagrangian evolution equations (2) are easily generalized to self-gravitating flows: 

f = g . (6) 

Here, the acceleration g is equivalent to the gravitational field strength according to 
Einstein's equivalence principle of inertial and gravitational mass. If we restrict ourselves 
to Newtonian theory, which is adequate for the considered range of scales, then g is 
constrained by the field equations: Vx X g = and Vx • g = A — A'kGq, where the 
density is again given by the deformation tensor similar to (4); for single-stream systems 
we have: Q = QqJ~^\ J '■= det[(9/i/9Xfc)]. We have to transform these field equations 
with respect to Lagrangian coordinates: 

9{gi], fp, fq) n ■ I'T J, ^ 

€pq[j—————— = ,t^J, {7a,b,c) 



d{Xi,X2,Xs) 
diga,fb,fc 



^ 2 d{Xi,X2,Xs) 

a,b,c 



J-^ = K-A-kGq . {Id) 



Inserting (6) and the density integral into (7) we obtain a set of four non-linear partial 
differential equations for the trajectory field f . Obviously, the gravitational problem is far 
more complex than inertial motion, which is determined by linear ordinary differential 
equations for f (eqs.(2)). (Alternative forms of eqs. (7) are to be found in Ehlers Sz 
Buchert 1994.) 

Currently, we are not able to give an exact solution to the equations (7) for generic 
initial data; we can pursue a perturbative analysis and solve the equations (7) up to a 
given order r with the following perturbation ansatz: 

r 

fM = a(0X + 5] sV^X,*) . (8) 

1=1 

The first term in this expansion represents the homogeneous and isotropic solutions of 
the system (7) (Friedmann-Lemaitre cosmogonies), in the simplest case a(t) oc t^^^ . 

Figure 3: The three— dimensional motion is illustrated in terms of an iso— density contour (with 
low density) at an instant of time, when most of the matter has moved from the sheet pattern towards 
filaments, and from the filaments towards the clump in the center of the picture (a £^4- singularity 
in the language of the 'Lagrange singularity theory'). In the 'pancake picture' this clump would be 
identified with the Coma cluster of galaxies embedded into a structure with the morphology of the 
'Great Wall'. 



III. 2. Relation to the inertial continuum 



Going to the first order in the perturbation approach (8) we obtain the following field of 
trajectories (Buchert 1992): 

f[i] = a{t) X + bi{t) U^(X)io + b2{t) U^(X)io + bs{t) W(X)i^ , (9) 

where U^, are the irrotational and rotational parts, respectively, of the initial 
peculiar-velocity (obtained by subtracting the initial Hubble-velocity of the homogeneous 
solution from the total velocity), and similarly W the peculiar-acceleration. A physically 
interesting class of first-order solutions is provided by the restriction to irrotational flows 
(U^ = 0), and to parallelity of velocity and acceleration (U^ = Wio)- This class is 
singled out by the growing mode solution of the Eulerian linear theory of gravitational 
instability. With this restriction we recover the well-known Zel'dovich approximation 
(Zel'dovich 1970, 1973); Zel'dovich suggested this extrapolation of the Eulerian linear 
theory into the weakly non-linear regime: 

f(^) =a(t) X + 6(t) U^(X)to . (10) 

The striking fact here is that we understand (10) as a perturbative solution of the 
full gravitational system (7), and immediately recognize a simple relationship to the 
inertial continuum studied thus far: the transformation to comoving deformations 
q = F(X,t) = a(t)~^f(X.,t) together with the transformation of the time variable 
t b{t)/a(t) plus the replacement of V(X) by U^(X) shows the formal equivalence 
of (10) and (3). All statements concerning generic structure elements which are valid in 
the case of inertial motion also hold to first order in the gravitational problem for this 
set of initial data. 

III. 3. Internal hierarchy of structures 

The perturbation approach (8) enables us to study higher-order gravitational effects of 
the structure formation process. At second order we observe a second shell-crossing 
singularity which arises as a consequence of the enhanced gravitational self-action of 
the three-stream system compared to the single-stream system: a second generation 
of pancakes is created inside the first pancakes (compare Buchert & Ehlers (1993) and 
my last Proceedings contribution, Buchert (1993), as well as Buchert (1994a) and ref. 
therein). Numerical N-body simulations show that this is only the onset of a hierarchy of 
shell-crossings which invoke a cascade of nested pancakes (see: Doroshkevich et al. 1980). 
We already noted that there exist three generations of pancakes due to the degeneracy 
of three eigenvalues of the deformation tensor in the case of inertial motion. 

Self-gravity creates a much richer differentiation of structures, which I call internal 
hierarchy of structures. In a one-dimensional submanifold this is illustrated in Fig. 4. 
Whether this 'self-similar', or 'fractal' substructuring (as suggested by Mandelbrot 1976) 
continues into a chaotic regime is unclear. There are indications that the first singularity 
already contains this infinity of bifurcations as is suggested by the study of the second- 
order solution (see the bifurcation diagram in Buchert Sz Ehlers 1993), and is supported 



by the behavior of the gravitational field strength near the singularity: a wide class of 
solutions to (7) admits 'first integrals' which have the structure that |g| oc g, i.e., the 
field strength blows up at caustics. This fact should be compared with the finiteness of 
g across caustics as found in Lagrangian perturbation solutions as well as in N-body 
simulations (see: Buchert 1994b). 




Figure 4: Same as Fig.l, however, in the case of self— gravitating motion, several foldings of jC, 
develop giving rise to several shell— crossings: a N + 2— stream— system (where N is the number of 
foldings) is enclosed by a hierarchy of nested caustics in Eulerian space; N is likely to approach OO. 
Note that here a one— dimensional spatial subsection is presented; the pure one— dimensional case, for 
which the first— order solution is exact (Buchert 1989), shows a single folding only. (The additional 
foldings have been added to Fig.l, therefore the single folding is still visible.) 



III. 4. External hierarchy of structures 



Up to now we studied the structure formation problem for initially coherent fields. In- 
deed, this picture can be used to successfully describe, at least qualitatively, present day 
observational data of the galaxy distribution: Identifying the central peak in Fig. 3 with, 
e.g., the Coma cluster of galaxies, then the filamentary connections between clusters and 
even sheet patterns can be found in the observations (the so-called Great Wall, Geller 
&L Huchra 1989). However, since sheet patterns have comparatively low densities, one 
usually only observes filamentary and clumpy patterns in the galaxy distribution, the 
sheets being hidden in the unseen dark matter constituent. 

There is considerable disagreement in the community as to whether this simplified 
picture is applicable; in other words, it is doubtful that such a coherence scale is present 
in the initial data at the recombination epoch. However, a pure baryonic universe as well 
as a 'Hot-Dark-Matter' universe show such coherence properties as a result of viscous, 
or free-streaming damping mechanisms, respectively. 

t 
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Figure 5l A scetch of the structure formation process in hierarchical models: Different spatial 
scales are marked; the collapse process is symbolized by crossing trajectories. Shell— crossing occurs 
earlier on smaller spatial scales. However, the merging of collapsed objects is controlled by the hier- 
archically superordered scale in such a way that the collapsed objects are "guided" by the collapse of 
the larger wavelength. Thus, truncation of higher frequencies in the initial power spectrum will not 
influence the structure formation process significantly on larger scales, as was demonstrated by Melott 
et al. 1994a,b. 

In order to study non-coherent initial data, we choose as description their power 



spectrum V{k), i.e., the ensemble average of the squared absolute values of the Fourier 
coefficients of the density field. For simplicity, we shall discuss cases in which the power 
spectrum has the form of a powerlaw: V= . The index n = —3 here plays a special 
role; it is that index for which the total power in a given fc— interval, J d^k 'P{k), is 
approximately constant (in fact logarithmic) as a function of scale k. 

In such a model, the self-similar linear growth of fluctuations implies that all scales 
enter the non-linear regime at about the same time. Consequently, all the structures 
we described in the present talk form on all scales simultaneously. For less small-scale 
power, we recover those models which have a sufficiently steep fall-off in the spectrum 
at short wave-lengths, thus forming structures on a scale with most power (usually the 
maximum of the spectrum like in a 'Hot-Dark-Matter' cosmogony). In the opposite 
case, we encounter those models which commonly are called hierarchical, i.e., structures 
form first on smaller scales and then merge into superunits on larger scales (as an ex- 
ample: the 'Cold-Dark-Matter' cosmogony). Especially the latter models result in that 
scale-dependent phenomenology which I call external hierarchy of structures (scetched 
in Fig. 5). Such models, in particular those which are close to n = —3 on small scales, 
are the most likely ones in comparison with observational data. That the picture of 
generic structure metamorphoses I described in this talk is applicable also in the case of 
hierarchical cosmogonies has been demonstrated by Melott et al. (1994a,b); see also the 
discussion in Buchert et al. (1994). 

IV. Open problems 

Summarizing, I have presented a picture which is characterized by structure metamor- 
phoses on all spatial scales with a continuing hierarchical feature of the internal structure 
differentiation down to smaller and smaller scales. 

This picture is limited on large scales by the horizon scale, where general relativistic 
models for self-gravitating motion have to be considered; on small scales, which are not 
ruled by the gravitational interaction, hydrodynamical efi'ects of the baryonic matter 
and chaotic thermal velocities of collisionless matter are relevant: still, self-gravitating 
structures idealized by images of generic singularities are thought to form the skeleton of 
large-scale structure down to the galaxy as structural unit. The gravitational interaction 
of multi-stream systems is an open problem in the description of structure formation 
behind caustics: the basic equations (7) have to be extended to the Vlasov-Poisson 
system in order to properly meet this problem. 

Another unsolved problem is the classification of generic structure elements in the 
case of self-gravitating motion: There is no proof specifying that class of motions which 
develop Lagrangian singularities. For general motions the mathematical background has 
to be extended. 

Since we observe hierarchical features on all spatial scales, it is not at all intruding 
that the universal structure admits a scaie of homogeneity, i.e., that the power spectrum 
bends over to a falling slope on the largest scales. It is likely that this hierarchy, which 
we observe in a large window of spatial scales, continues also to scales which are be- 
yond the access of current observations. Indications from the COBE measurement of the 



large-scale end of the power spectrum that this is not the case I consider as preliminary: 
cosmic variance arguments express the statistically incomplete status of this observation 
(see: Scaramella & Vittorio 1993). Both internal and external hierarchies may be an 
indication that we have to develop a fully hierarchical model of the universe abandon- 
ning the prejudice that the universe should follow the Priedmann-Lemaitre cosmogonies 
on average and that it should admit scale-independent average characteristics (like the 
Hubble constant) calculated from these simple-minded models. 

Acknowledgements: 

I would like to thank Matthias Bartelmann for providing Fig.3 as well as Thomas Boiler 
for scanning the figures. I am greatful to Profs. Xiaoyang Xia and Zugan Deng for their 
efforts to make this workshop most enjoyable. 

References 

Arnol'd V.I. (1982): Trudy Sem. Petrovskogo 8, 21 (in Russian). 

Arnol'd V.I., Shandarin S.F., Zel'dovich Ya.B. (1982): Geophys. Astrophys. Fluid Dyn. 
20, 111. 

Arnol'd V.I., Gusein-Sade S.M., Varchenko A.N. (1985,1988): Singularities of Differen- 

tiable Maps Vol. I, II, Birkhauser, Boston. 
Buchert T. (1989): Astron. Astrophys. 223, 9. 
Buchert T. (1992): M.N.R.A.S. 254, 729. 

Buchert T. (1993): Proc. 4. MPG CAS workshop, Ringberg, Tegernsee; eds.: G. Borner, 

T. Buchert, MPA/P8, 204. 
Buchert T., Ehlers J. (1993): M.N.R.A.S. 264, 375. 
Buchert T. (1994a): M.N.R.A.S. 267, 811. 
Buchert T. (1994b): Phys. Rep., submitted. 
Buchert T., Shandarin S.F., Weifi A.G. (1994): in preparation. 

Cao S.-L. (1993): Proc. 4. MPG CAS workshop, Ringberg, Tegernsee; eds.: G. Borner, 

T. Buchert, MPA/P8, 193. 
Doroshkevich A.G., Kotok E.V., Novikov I.D., Polyudov A.N., Shandarin S.F., Sigov 

Yu.S., (1980): M.N.R.A.S. 192, 321. 
Ehlers J., Buchert T. (1994): M.N.R.A.S., to be submitted. 
Geller M.J., Huchra J.P. (1989): Science 246, 887. 

Mandelbrot B. (1976): Comptes Rendus (Paris) 282A, 119 (in French). 
Melott A.L., Pellman T.F., Shandarin S.F. (1994a): M.N.R.A.S. 269, 626. 
Melott A.L., Buchert T., Weifi A.G. (1994b): Astron. Astrophys., in press. 
Scaramella R., Vittorio N. (1993): M.N.R.A.S. 263, L17. 
Shandarin S.F., Zel'dovich Ya.B. (1989): Rev. Mod. Phys. 61, 185. 
Thom R. (1975): Structural Stability and Morphogenesis, W.A. Benjamin Inc., Mas- 
sachusetts. 

Zel'dovich Ya.B. (1970): Astron. Astrophys. 5, 84. 
Zel'dovich Ya.B. (1973): Astrophysics 6, 164. 

Zel'dovich Ya.B., Myshkis A.G. (1973): Elements of Mathematical Physics - Medium of 
Non-interacting Particles, Nauka, Moscow (in Russian). 



